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Abstract 

We study the early-stage kinetics of thermodynamically unstable 
systems with quenched disorder. We show analytically that the growth 
of initial fluctuations is amplified by the presence of disorder. This is 
confirmed by numerical simulations of morphological phase separation 
(MPS) in thin liquid films and spinodal decomposition (SD) in binary 
mixtures. We also discuss the experimental implications of our results. 
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Introduction 



There has been intense research interest in the kinetics of phase transitions in 
systems which have been rendered thermodynamically unstable by a sudden 
change of parameters, e.g., temperature, pressure, magnetic field, etc. The 
subsequent far-from-equilibrium evolution of the system is usually nonlinear, 
and is characterized by complex spatio-temporal pattern formation. The sys- 
tem approaches its new equilibrium state via the emergence and growth of 
domains enriched in the preferred phases. This nonequilibrium dynamics is 
often referred to as phase- separation kinetics or domain growth or coarsening 
or phase ordering dynamics [TJ. These processes occur over a wide range 
of length-scales and time-scales, and are of great scientific and technological 
importance. 

In pure (disorder-free) systems, the coarsening domains are usually char- 
acterized by a divergent length scale which grows as a power-law, L(t) ~ 
where the exponent <fi depends on the transport mechanism. For example, in 
the phase separation of an unstable binary mixture via diffusion, L(t) ~ t 1//3 
in the late stages of evolution. This is known as the Lifshitz-Slyozov (LS) 
growth law. In recent works J2j[3], we have established that the LS growth law 
also describes the late-stage dynamics in the morphological phase separation 
(MPS) of an unstable thin liquid film (< 100 nm) on a solid substrate. During 
MPS, the film segregates into flat regions (domains) and high-curvature re- 
gions (defects or hills). In these works, we have also emphasized the analogies 
and differences between the kinetics of phase separation in binary mixtures 
and thin films. 

Of course, real experimental systems are never pure or disorder-free: 
they are invariably characterized by chemical and physical heterogeneity 
HI El El [71 [H El HUl HH H21 H3]. Therefore, it is natural to investigate the 
effects of disorder on the domain growth process. For segregating AB mix- 
tures, coarsening is driven by interfaces between coexisting v4-rich and 5-rich 
domains. These interfaces are trapped by sites of quenched disorder, and 
subsequent growth occurs by thermally-activated hopping of interfaces over 
disorder traps [H]. Thus, domain growth is drastically slowed down in the 
late stages, though the domain morphology is unaffected [151 HSJ [TO EE]- On 
the other hand, unstable liquid films undergoing MPS segregate into regions 
with heights h = h m and h = oo, i.e., there are no coexisting phases. The 
"interface" between the flat phase {h = h m ) and the high-curvature steepen- 
ing phase {h = oo) becomes sharper as time proceeds, and is unaffected by 
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the presence of disorder. Thus, the asymptotic regime of MPS with disorder 
is universal, and is characterized by the LS growth law [3J. 

In this letter, we focus on the effect of disorder on the early stages of phase 
separation, i.e., the growth of small fluctuations about a homogeneous initial 
state. As we discuss shortly, both binary mixtures and thin films [19] are de- 
scribed by the Cahn-Hilliard-Cook (CHC) equation [201 EI] with appropriate 
free energies. The early-stage dynamics is described within the framework 
of a linear theory, usually referred to as CHC theory. In this letter, we re- 
port that growth in the CHC regime is strongly amplified by the presence of 
disorder. We present both analytical and numerical results to support this 
scenario. This should be contrasted with the late-stage dynamics, which is 
either slowed down (for binary mixtures) or unaffected (for thin films) by the 
presence of quenched disorder. 



Dynamical Equations and Analytical Results 



The starting point of our discussion is the dimensionless form of the CHC 
equation, which describes the evolution of a conserved order parameter tp(f, t) 
at space-point f and time t [22] : 



0_ 

dt 



ip(r,t) = V ■ 



M(V>)V 



SF\ 



Sip 



+ rf(r, t) 



(1) 



In Eq. (CQ), M(i/j) is the mobility, and the Gaussian white noise ff satisfies the 
fluctuation-dissipation relation: 



fftr, t) 



7fc(f',t%(r",t") 



0, 

2eM(ip)5 ij 5(r' - r")S(t' -t"). 



(2) 



Here, the bars denote an averaging over the noise ensemble, and e measures 
the strength of the noise. Further, F[i/j] is the free-energy functional which 
has the following form: 



dr 



(3) 



In Eq. ()3]), fa(i>) denotes the local free energy which is parameterized by the 
disorder strength a. The second term on the RHS refers to the interfacial 
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energy or surface tension. Replacing Eq. (E]) in Eq. flTJ, we obtain 

^(f,t) = V-[M(^)V{f'M-^}+v(r,t)] . (4) 

The dimensionless Eqs. dX|)- (J4j) are obtained from the corresponding di- 
mensional equations by using the natural scales of length, time, and order 
parameter pQ. Before proceeding it is useful to discuss the functional form 
of f a {ip) for the problems of interest here. For unstable thin films, ip is the 
height (usually denoted as H), and a typical form of /o(V0 f° r the pure case 
(a = 0) is 



fo(H) = -i 



1 - R R 

+ "T7n 



l(H + D) 2 ' # 2 ■ (5) 

This potential describes a thin film on a coated substrate. In Eq. ([5]), is 
the ratio of the effective Hamaker constants for the system, and D is the di- 
mensionless coating thickness. One can introduce the disorder in the coated 
potential through R (chemically heterogeneous substrate) or D (physically 
heterogeneous substrate). Here, we consider the case with chemical hetero- 
geneity, so that R is a r-dependent random variable uniformly distributed in 
the interval [R m — a, R m + a}. 

For segregating binary mixtures, ip denotes the density difference of the 
two species: 0(r, t) = n^(r, t) — rie(r, t), where n a is the density of species a. 
The local free energy for the disorder-free case is modeled by a 4 -potential 
with a double- well structure: 

+ £ >=i- (6) 

We model chemical heterogeneity by making b a space-dependent random 
variable, which is uniformly distributed in [1 — a, 1 + a], i.e., the critical 
temperature varies from point to point. 

Let us now return to Eq. (J3J). To understand the dynamics of the early 
stages, we consider the growth of initial fluctuations, ip(r, t) = ip + S'p(f, t). 
Here, ipo is the homogeneous value of the order parameter in the unstable 
state. The linearization of Eq. 03]) yields 

J^(r,0 = M(^ )V 2 (-a - V 2 )^ + V • ff(f,t), (7) 

where we have introduced the parameter a = —f"(i>o)- We neglect terms 
involving the spatial derivatives of the disorder variable, as these are only 
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relevant on microscopic scales. (We will subsequently average all physical 
quantities over the disorder distribution.) The Fourier transformation of 
Eq. © (with wave-vector k) gives 

d -> 

— 5<4>(k, t) = M(ij Q )k 2 (a - k 2 )5ij(k, t)-ik- ff(k, t), (8) 

which has the following solution: 

5ip(k, t) = exp [M(ij )k 2 (a - k 2 )t] 5ip(k, 0) 

-ik- [ dt'exp[M(^ )k 2 {a-k 2 )(t-t')]f]{k,t'). (9) 
Jo 

To obtain Eq. we have set M(i[)) ~ M(ip ) in Eq. ([2]), which is valid at 
early times. Clearly, averaging over the noise ensemble yields 

ty(k, t) = exp [M(ij )k 2 (a - k 2 )t] 5ip(k, 0). (10) 

The parameter a is spatially uniform for a chemically homogeneous sys- 
tem. For the disordered system, we have the local free energy 

/«(V0 = /o(V0 + z/i(V0, (ii) 

where x is a uniformly-distributed random variable in the interval [—a, a]. 
For the disorder types we consider here, Eq. (Hip is exact. Otherwise, it 
represents the first two terms of a Taylor expansion in the disorder x. Thus, 

Finally, integrating Eq. f JTOj) over the uniform disorder distribution with 
the probability P(x) = l/(2a) gives 

ty(k,t) = J dx^-exp [M(i) )k 2 {a(x) - k 2 }t] 5ip(k,0). (12) 

Some simple algebra results in the following expression for the disorder- 
averaged growth of fluctuations: 

Tilt +\ rMY/M.2/ f/ir i \ 1 2\ + i sinh [ M ('4 ) o)fi(ipo)k 2 at} r 

5^k,t)=e W [M(^)k {-f M-k }t] mM{ ^ )k2at *KM). 

(13) 

The factor sinh(ct) / (ct) determines the effect of disorder on the order-parameter 
evolution. This factor tends to 1 as a — > or k — > or t — > 0. However, 
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the growth of initial fluctuations can be strongly amplified for nonzero a as 
sinh(ct) grows exponentially with t. 

Finally, we evaluate the time- dependent structure factor, which is probed 
in scattering experiments using, e.g., light, neutrons. This is defined as 
follows: 

S(k,t) = {6ip{k,t)6ip*{k,t)), (14) 

where the angular brackets indicate an averaging over thermal noise and 
independent initial conditions. Using the expression for Sijj(k, t) from Eq. (Q, 
we obtain 

S(k, t)= AV exp [2M(tfj )k 2 (a - k 2 ) t] 

+ T^rp) i ex P i 2M (^)k 2 (a - k 2 ) t] - 1} . (15) 

To obtain Eq. (fl5|) . we assume that {S^(r, 0)5ip*(r', 0)) = A5{f — r'), where 
A is the amplitude of initial fluctuations in the order parameter. Further, 
V denotes the volume of the system. The corresponding disorder-averaged 
structure factor is 

§(M) = A Vexp [2MW)e{-XW) ~ ^}t] "t^'t^'tw" 

2M{il) )f'{(i) )k 1 at 

V i-fUm-afl' (^o)-fc 2 i exp [2M(iP )k 2 yt] - 1 1 



Numerical Results 

Next, we describe our numerical results for the chosen systems. For MPS in 
a thin film, the appropriate model is Eq. (j3J) with ip = H and the local free 
energy: 



f a (H) = f (H)+xf 1 (H), h{H) = \ 

o 



1 



[H + D) 2 H 2 



(17) 



where fo{H) is defined in Eq. (JSJ). Here, we consider the case with mobility 
M(H) = H 3 (corresponding to Stokes flow with no slip) and e = 0, i.e., 
without noise. We numerically solve Eq. (j3J) in d = 1 starting with a small- 
amplitude (~ 0.01) random perturbation about the mean film thickness, 
H = 1. The system size is uLm, where Lm is the dominant wavelength 
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for R = R m (n ranges from 16 to several thousands). We apply periodic 
boundary conditions at the end points. A 512-point grid per L M was found 
to be sufficient when central-differencing in space with half-node interpolation 
was combined with Gear's algorithm for time-marching, which is convenient 
for stiff equations. The parameters D = 0.5, R m = —0.1 and a were chosen 
so that the film is spinodally unstable at H — 1 (i.e., a > 0) for all values of 



For segregation in a disordered binary mixture, we consider Eq. (j3J) with 

if) = (f) and the free energy: 



where fo(<f>) is defined in Eq. (j6]). We consider the case with M(</>) = 1 and 
e = 0. We numerically solve Eq. (j3J) via an Euler-discretization scheme on a 
d = 2 lattice of size 512 2 . The mesh sizes in space and time are Ax = 0.5 and 
At = 0.001, respectively. The spatial mesh size is small enough to resolve the 
interface region between A-rich (0 = +1) and S-rich (0 = —1) domains. The 
interface thickness £ ~ a/2 in our dimensionless units. The initial condition 
for a run consists of small- amplitude (~ 0.01) fluctuations about 0o — 0) 
corresponding to a critical AB mixture with 50%A and 50%5. 

First, we present results for MPS in films with chemical heterogeneity. 
Recall that we are interested in the effect of disorder on the early stages of 
growth, i.e., the CHC regime. The height profiles at a fixed time (t = 4400) 
for various disorder strengths are shown in Fig. [T^a). The initially random 
perturbations grow exponentially in time, with fastest growth for the case 
with maximum disorder (a = 0.04). The growing fluctuations rapidly select 
the most unstable wavelength, A = Lm - the location of the zero crossings 
is unaffected by the disorder. However, the peaks of the profiles are strongly 
amplified due to the factor smh(ct) / (ct) in Eq. (fl3|) . This amplification may 
be quantitatively characterized by studying the statistical properties of the 
height profile. In Fig. CUb), we plot C(0, t) = (5H(x,t) 2 ) vs. t. This quantity 
is related to the structure factor in Eq. ( JT4l) as 



R. 



/a(0)=/ O (0)+xA(0), 




(18) 




(19) 



where d is the dimensionality. As expected, C(0, t) grows fastest for the case 
with maximum disorder. 
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The initial exponential growth of the profiles in the CHC regime is fol- 
lowed by a nonlinear saturation which results in the formation of domains of 
the flat phase, separated by defects of the high-curvature phase. These do- 
mains grow via the diffusive transport of material. The CHC regime and the 
domain growth regime are both shown in Fig. [2j where we plot the number 
of defects (hills) N(t) vs. t. The characteristic length scale is L(t) ~ iV(£) _1 . 
The CHC regime is followed by an intermediate stage where disorder sites 
trap the emergent interfaces. In this regime, the more disordered systems 
show slower growth. The intermediate stage is very sensitive to the presence 
of even small amounts of disorder. In the late stage, the domain growth is 
unaffected by disorder. In this regime, the interfaces (between H = H m and 
H = oo) have steepened too much to be captured by the local disorder [3J. 
Thus, this regime shows universal LS growth [N(t) ~ i -1 / 3 or L(t) ~ t 1 / 3 ]. 
This result has important experimental implications, viz., it is not necessary 
to work with specially-prepared extra-pure substrates, at least in the context 
of domain growth studies. 

Next, we present our results for the early stages of spinodal decomposition 
(SD) in d — 2 binary mixtures. In Fig. [3J, we show the variation of the order 
parameter along a diagonal of the lattice (y = x). We plot profiles for the 
pure (a = 0) and disordered (a = 0.5, 0.8) cases at t = 20. Again, we see that 
there is an amplification of fluctuations for higher disorder values, though the 
selected wavelength does not change. 

Finally, we study the structure factor S(k, t) which provides a quantitative 
measure of the phase-separating morphology. In Fig. 0J we plot S(k, t) vs. 
k at t = 20 for the pure and disordered cases shown in Fig. |3j In the 
early stages or the CHC regime, the structure factor grows exponentially 
in time but the peak position does not shift, i.e., there is no change in the 
characteristic scale. We see that the fastest growth occurs for the highest 
disorder amplitude. The solid lines in Fig. H] denote the expression in Eq. ( TT6|) 
with e = 0. Our numerical results are seen to be in excellent agreement with 
theory. 

In the late stages of phase separation, disorder has the opposite effect, 
drastically slowing down domain growth [T5J [TBI CEO IB] • As stated earlier, 
this is because interfaces between coexisting A-rich and 5-rich domains are 
trapped by sites of quenched disorder. This should be contrasted with the 
universal late-stage scenario for segregating films, as described above. 
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Summary and Discussion 



In summary, we have studied the effect of quenched disorder on the early 
stages of evolution in two important systems: unstable thin films and phase- 
separating mixtures. Both these systems are described by the Cahn-Hilliard- 
Cook (CHC) equation with appropriate choices of the free energy. We find 
that disorder amplifies the exponential growth of fluctuations in the early 
CHC regime in both systems, and we have obtained a simple analytical ex- 
pression for the amplification factor. In the asymptotic regime, there is a 
major difference between disordered domain growth in these two systems, 
viz., films undergoing MPS are unaffected by disorder, whereas segregating 
mixtures are drastically slowed down by disorder. Our results have impor- 
tant experimental consequences, and we hope that they will be subjected to 
experimental tests. 
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Zero disorder 




Figure 1: Early stages of morphological phase separation (MPS) in an un- 
stable thin film with quenched disorder, (a) Height profiles [H(x, t) vs. x] at 
t = 4400 for different disorder amplitudes. The simulation details are pro- 
vided in the text, (b) Plot of (5H(x,t) 2 ) vs. t for the height profiles shown 
in (a). 
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Figure 2: Variation of the number density of defects with time [N(t) vs. t] 
for MPS in an unstable thin film. Curves 1 to 6 present results for disor- 
der amplitudes a = 0, 0.0001, 0.0004, 0.002, 0.01 and 0.04, respectively. The 
dashed lines with slopes of —1/4 and —1/3 denote exponents for the early 
and late stages, respectively [2j [3] . 
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Figure 3: Early stages of spinodal decomposition (SD) in an unstable disor- 
dered binary mixture in d — 2. The plot shows the variation of the order 
parameter along the diagonal [<f>(x, x, t) vs. x] at t — 20 for different disorder 
amplitudes. The simulation details are given in the text. 
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Figure 4: Plot of the spherically-averaged structure factor [S(k,t) vs. k] at 
t = 20 for different disorder amplitudes. The structure factor is obtained 
as an average over 50 independent runs. The initial condition for each run 
consists of a different order parameter and disorder configuration. The solid 
lines denote the expression in Eq. (fT6|) with e = 0. 
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